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ABSTRACT 



Basic heat transfer empirical and analytic equations are 
applied to a double envelope airship concept which uses 
heated Helium in the inner envelope to augment and/or 
control gross lift. The convective and conductive terms 
lead to a linear system of five equations for the concept 
airship, with the nonlinear radiation terms included by an 
iteritive solution process. 

A FORTRAN program is used to perform the tedious 
calculations, and graphed results are obtained for the 
variables of interest. These indicate that a simple use of 
airship engine exhaust heat will give more than a 30% 
increase in gross airship lift. Possibly more than 100% 
increase can be achieved if a "steam injection" heating 
system, with associated design problems, is used. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 



PROGRAM 

NAME 


TEXT 

SYMBOL 


DEFINITION 


A 


A 


Generic area or matrix value 


AHTR 


^htr 


Heater effective radiant area 


AIN 




Inner envelope effective area 


AOUT 


^out=^o 


Outer envelope area 


B 


B 


Generic matrix value 


C 


C 


Generic matrix value 




''Pair 


Specific heat of air 




d 


Envelope spacing 


DELT 


Delta T 


Temperature difference 


DLH20 




Change in lift due to steam volume 


DLIFT 




Change in lift (Ibf) from temperature 
increase 


DLTOT 




Total change in lift 


DMDT 




Equivalent fuel burn rate 


DQDTN 




"New” value of heat flux computed in 
present iteration 


DQDTRO 


^rO-2 


Radiative heat flux from heater surface 
to the inner envelope 


DQDTR3 


<3r3-4 


Radiative heat flux between the 
envelopes 


DVH20 




Increase in volume due to steam 


E 


E 


Generic matrix value 


EHTR 


^htr 


Emissivity of heater surface 


EIN 


£-in 


Emissivity of inner envelope 


EOUT 


£ out 


Emissivity of outer envelope 
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GR 


Gr 


Grashof Number 


GRD 




Grashof Number based on envelope 
radius difference 


GRFL 


Grfi 


Grashof Number based on fluid 
temperature 




h 


Convective heat transfer coefficient 


HI 


hi 


Convective heat transfer coefficient 
for the inner volume 


HIAIR 


hair 


Convective heat transfer coefficient 
based on properties of air 


H4 


h4 


Convective heat transfer coefficient 
for the outer surface 


KAIR 


^air 


Thermal conductivity for air 


KEQV 


’^eqv 


Thermal conductivity- 
equivalent computed value 


KFL 


^fl 


Thermal conductivity 
based on fluid temperature 


KHE 


(D 


Thermal conductivity for Helium 


KIN 


hin=^^i 


Thermal conductivity 
for the inner envelope 


KINF 


^inf 


Thermal conductivity (of air) 
based on outside air temperature 


KMEAN 


hmean 


Thermal conductivity (of air) 
based on volume mean temperature 
between the envelopes 


KOUT 


^out~^o 


Thermal conductivity 
for the outer envelope 


L 


L 


Airship length 


MUAIR 


^air 


Viscosity of air 


MUFL 


Afl 


Viscosity based on fluid temperature 


MUHE 


(D 


Viscosity of Helium 


MUI 


Ainf 


Viscosity based on outside temperature 
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MUMN 




Viscosity based on volume mean air 
temperature 




Nu 


Nusselt Number 




N^air 


Nusselt Number for air 






Nusselt Number for Helium 


PDL 




Percent change in lift due to 
temperature increase only 


PDLTOT 




Total percent change in lift 




Pr 


Prandtl Number 


PRFL 


Prfi 


Prandtl Number based on fluid 
temperature 


PR I 


P^inf 


Prandtl Number based on outside air 
temperature 


PRMN 


P^mn 


Prandtl Number based on mean fluid 
temperature 


Q 


q 


Heat flux 


RE 


Re 


Reynolds number 


REL 


(1) 


Reynolds number based on airship length 


RHOAIR 


P air 


Density of air 


RHOFL 


Pfl 


Density based on fluid temperature 


RHOHE 


pHe 


Density of Helium 


RHOI 


P inf 


Density based on outside air 
temperature 


RHOMN 


Pmn 


Density based on mean fluid 
temperature 


RIN 


Pin 


Inner envelope effective radius 


RINIT 




Iteration initial trial radius 


RMAX 


Pitiax 


Maximum radius 


RMEAN 


Rm 


Average of inner and outer envelope 



radii 
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ROUT 


^out 


Outer envelope radius 


RTRY 




Present iteration trial radius 




shp 


Shaft horsepower 


SPV 


V 


Specific volume (of water vapor) 


SUM 




Value used in system solution 


T 


T 


Generic temperature 


TO 


To 


Inner volume mean fluid temperature 


T2 


T2 


Temperature on inner surface of 
inner envelope 


T3 


T3 


Temperature on outer surface of 
inner envelope 


T4 


T4 


Temperature on inner surface of 
outer envelope 


T5 


T5 


Temperature on outer surface of 
outer envelope 


TFL 


Tfl 


Fluid temperature 


THE 


<D 


Helium temperature 


THIN 




Thickness of inner envelope 


THOUT 


to 


Thickness of outer envelope 


THTR 


Thtr 


Heater surface temperature 


TI 


Tinf 


Temperature of outside air 


TMEAN 


Tm 


Volume mean temperature of the air 
between the envelopes 


UI 


Uinf 


Velocity of the airship 


VH20H 




Volume of the water vapor with the 
Helium heated 


VH20I 




Volume of water vapor with the Helium 
unheated 


VIN 


^in 


Inner envelope volume 


VTOT 


^out 


Outer envelope, or total volume 
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I. 



NATURE OF THE PROBLEM 



A. CONTROL OF LIFT 

In a fixed or rotary wing aircraft, the control of lift 
is normally a simple problem of changing engine power and/or 
angle of attack on the wing or rotor, and desired lift is 
quickly attained up to the aircraft's capability limits. 
Lighter than air ships, however, have a more inherent 
problem in increasing or decreasing lift. Typically, 
changes are produced by adding or valving lifting gas, 
ballasting by some heavy and readily obtained material such 
as sand or water, or employing engine exhaust based water 
recovery systems to compensate for fuel burnoff. 
Disadvantages of these are well known: availability and 

cost of lifting gas; ballasting material is not always 
available or easily transferred, and is a one way operation 
once airborne; and water recovery systems are bulky and 
heavy. 

Control of the lifting gas temperature could provide a 
more flexible and complete means of lift control. In 
addition, a significant increase in lifting gas temperature 
will ideally provide a proportionate increase in gross lift. 

B. HISTORY 

Using heated lifting gas to control or augment the gross 
lift of an airship is an old idea, and an obvious extension 
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of early hot air ballooning. Burgess, in his book "Airship 
Design" , lists this idea in his Airship Fallacies chapter 
due to his estimation on the rate of heat loss. He does, 
however, offer "...gas cell construction in which convective 
losses were greatly reduced might possibly result in a 
practical solution of the problem" (reference 1, pp. 288- 
289). Present day fabrics, films and laminates offer much 
more flexibility in such construction than the rubberized 
cotton in use in the early 1900's (references 1 and 2). 



C. RESEARCH FINDINGS 

A library and information search was conducted to 
determine whether research results were available for such a 
heat transfer problem. The following results were obtained: 

- Mention in NASA CR137692 (reference 3) that Airfloat 
Transport, Ltd., had looked at heated Helium at takeoff 
to eliminate water recovery apparatus. No further 
information was obtained. 

- Concepts on using engine heat to heat lifting gas, but 
no technical evaluation by Davenport (reference 4) . 

- A letter from Mr. I. P. Alexander of Wren Skyships, 

Ltd. , stated that Wren Skyships had looked at ballonet 
heating on the RS-1. Additional mention was made of 
electrical heating tested on the LZ-130 Graf Zeppelin 
II. Neither of these concepts proved feasible, but 
these involved a small modification, and not a basic 
design concept. No technical data was received. 

- In a telephone conversation with Mr. Mark West of 
Aerostar, he stated that Aerostar had worked with 
students at Rensselaer Polytechnic Institute on a steam 
lift blimp design. No technical results were 
obtainable. 

- Extensive hot air and high altitude balloon heat 
transfer modelling is available. Applicable data from 
Stefan (reference 5) was modified for use in this model. 
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Also, applicable gas and film radiative properties were 
investigated in reference 6. 

The conclusion based on this available information is 
that only limited applications of lifting gas heating have 
been seriously explored. If a design dedicated to the use of 
heat to control and augment the lift of a Helium airship has 
been completed, no comprehensive mathematical modelling for 
such a design is available. 
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II. THE PHYSICAL MODEL 



A. A PRACTICAL MODEL 

A comprehensive design concept for an airship is 
certainly beyond the scope of this paper. It is, however, 
desirable to have a reasonable physical concept on which to 
base the mathematical modelling and derivations. That 
concept, and some promising alternatives, are briefly 
explored, and the basis for the mathematical model is set 
forth. 

As encountered in any design, conflicting requirements 
must be considered, and the final choice of the model is an 
optimized compromise. Reduction of convective heat 
transfer, and maximum resistance to heat loss requires a 
multiplicity of panels within a volume, and insulative 
material coated onto panels. It is quickly realized, 
however, that such an airship will not be "lighter than 
air", and a minimum of envelope yardage and insulation must 
be pursued. 

In the absence of design specifications and "tradeoff" 
parameters, and in order to generate a simple but effective 
hypothetical physical model, it is assumed that the weight 
penalty of insulation material in addition to the inherent 
insulation properties of the envelope fabrics is 
prohibitive. Likewise, the use of additional fabric panels. 
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other than a basic double envelope design, is assumed 
prohibitive. The resulting concept is basically a 
conventionally shaped airship (see Figures 1 and 2) with the 
outer envelope enclosing lifting gas and ballonet volumes, 
providing basic protection from the external environment, 
aerodynamic shape, and pressure containment to hold that 
shape . 

An inner envelope contains the lifting gas, and must be 
expandable to accommodate both heating and altitude 
variations required for the airship. This envelope, then, 
must also carry the weight of the airship and loads. The 
weight suspension system also needs to constrain the shape 
of this inner envelope in such a way that the gap between 
inner and outer envelopes is reasonably uniform in order to 
minimize heat loss. The concept depicted in Figure 2 is of 
an adjustable/variable catenary cable system suspending the 
airship loads from the top of the inner envelope. The 
lengths of the cables are varied through a winch and pulley 
system to hold a desired pear shape, or natural buoyancy 
form, within the outer envelope. 

The air space between the envelopes is used as the 
ballonet in a standard fashion, and could, of course, be 
subdivided as required for safety and airship control . This 
volume also proves to be the most effective heat loss 
barrier as might have been expected. 
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Figure 1 : Airship Concept Geometry Overview 





Figure 2: Inner Envelope Suspension 
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In order to get a solid "feel" of the results from the 
heat loss rate calculations, specific external dimensions 
and gross lift capability were chosen to coincide closely 
with the probable sizes of near future NAVY or Coast Guard 
airships. These are: 

- gross lift of 130,000 Ibf (2 x 10^ ft^ of He) 

- length of 398.9 ft 

- maximum diameter 133 ft 

- slenderness ratio (L/2Rmax) = 3 

- total volume = 4 x 10^ ft^ 

Note that the outer envelope volume is twice the initial 
inner envelope volume. 

B. HEATING CONSIDERATIONS 

A practical, high capacity, efficient and lightweight 
system to provide heat into the inner envelope volume is 
critical to this concept. Low heat intensive systems such 
as heat pumps are assumed to be prohibitively heavy in the 
capacity ranges required. A simple air fed and vented 
petroleum based liquid burning heater may be a reasonable 
solution to these requirements. 

Davenport, in reference 4, proposes the use of engine 
exhaust to heat this volume. Appendix C shows a simple 
derivation of the amount of waste heat produced in engines 
such as may be used for propulsion of an airship this size. 
The problem of efficiently transferring this heat to the 
inner volume still remains. Forced air ducting and/or 
exhaust ducting may be the best, if least efficient means. 
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since the heat exchange system required between hot 
ducting or heater surfaces and the inner Helium volume might 
be extensive and heavy, the injection of steam produced from 
a low pressure boiler or engine manifolds directly into the 
Helium is considered as a possible alternative. Section 
of this paper contains some comments on the 
difficulties of estimating the inner envelope heat transfer 
characteristics with saturated Helium and the condensation 
of water vapor on the inner envelope. Some foreseen 
detrimental factors in the use of such a system are the 
weight of the condensate on the inside of the inner 
envelope, corrosion or environmental effects on the inner 
envelope, weight of water required, and the energy to 
evaporate that water. A rough calculation of these values 
is completed in Appendix D. Benefits of the use of such a 
system are seen in section IV. 

C. ALTERNATIVES AND IMPROVEMENTS 

Using the inner volume as ballonet, with the lifting gas 
volume between the envelopes, has two significant advantages 
over the chosen model. First, the more elaborate suspension 
system for the inner envelope is not required. Second, an 
open burner or direct injection and mixing of heated air 
into the inner volume makes a heat exchange ducting network 
unnecessary. A third benefit, heating the heavier air, and 
thus reducing the density proportionately, is seen to be a 
fallacy due to the fact that a fixed Helium mass is 
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retained, while the air mass is discharged to hold a 
constant total volume. This alternative concept has the 
disadvantage of a reduced heated volume as temperature or 
altitude are increased, while in the chosen concept, the 
directly heated volume will increase with these parameters. 

An easy improvement to the modelled concept would be the 
use of lightweight panels such as a mylar film to break up 
the convective patterns. No estimate of the effectiveness 
of these panels is attempted in the present model. 

Further research might result in fabrication techniques 
for envelopes with dramatically increased heat loss 
resistance. One approach is to trap relatively thick air 
layers between mylar films of the fabric composites, which 
could result in a factor of 10 increase in this value. 
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III. THE MATHEMATICAL MODEL 



A. THE HEAT TRANSFER PROBLEM IN GENERAL 

Analytical solutions to the convective heat transfer 
problem are limited to laminar flow over specified 
geometries (reference 7) . Empirical data curve fits must be 
used to apply known heat transfer characteristics to this 
airship heat loss problem. Since the available empirical 
equations apply to specific geometries, the airship shape is 
initially idealized as in Figure 3, with a hemispherical 
nose, cylindrical body, and conical tail. 

With a fixed slenderness ratio for the external shape, 
L/2Rmax = 3, the relationship between maximum diameter, 
envelope surface area and volume in this idealized model is 
fixed from basic solid geometry as: 

Rj^ax = (3V/13/7)1/^ = (A/ (8+/5) ) V2 

A= (8+i/5) (3V/13/7 ) 2/3 

The inner envelope is idealized as a shape concentric to 
the outer envelope (Figure 3) , with the total length the 
same as the outer envelope, but with the diameter varied as 
required to contain the lifting gas volume with changes in 
temperature (and altitude) . With a constant proportion for 
the hemispherical nose and conical tail, the required length 
is derived for the cylindrical body, and again from solid 
geometry : 
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Figure 3: Idealized Envelope Shape 



''in = ril. . 5T7RinV3 

Ain “ 2/7Rin(L - 3Rij^) 

In these idealized geometries, volume and area are 
greater than typical airship envelopes for a given maximum 
diameter. The more important volume and area relationship, 
however, is very close. These relationships are used 
because of their simplicity, allowing quick computer 
iterations, and the variation in manufactured airship 
envelope shapes and their corresponding geometric 
relationships . 

B. CONDUCTION AND CONVECTION MODELLING 
1 . Conduction Through the Envelopes 

Since the dimensions of the airship are very large 
compared to the thickness (t) of the envelopes, a one 
dimensional heat conduction equation is completely valid 
(reference 8, ch. 1): 

q = -(kA/t) (T^^t - Tin) 

The envelope thickness is idealized as a homogeneous 
material, although Vadala, in reference 2, clearly expects a 
multilayered fabric composite, such as a Kevlar base with 
covering Mylar films. The thermal conductivity (k) values 
used in this model are derived from reference 9 based on 
properties of Kevlar fabric in an appropriate thickness of 
.05 in. 

k = 0.03266 BTU ft/hr ft^ F 
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External Convection 



2 . 

Depending on operating conditions, the external flow 
will meet the criteria for forced, mixed, or free convection 
heat transfer conditions. In this approach, only forced and 
free convection values are derived, and a curve is sketched 
between zero and ten knot forced velocity values to 
represent the mixed convection regime. 

Operationally, this regime is expected to prove very 
important, as heat up and cool down periods will be required 
on takeoff and landing, depending on airship loading. With 
the airship tethered, light winds or even eddy currents in 
relatively calm air, with the boundary layer effects of the 
earth's surface, are assumed to make a direct analysis of 
this regime inaccurate. Therefore, even though a limited 
analysis of this mixed flow regime is possible, only the 
rough curvefit is used in the present model. 

For the free convection condition, the airship is 
idealized as a horizontal cylinder of equivalent area, and 
equation 7-36 from reference 8 is chosen as the closest 
empirical fit to the physical conditions: 

NUf = hL/k 

= {0. 6+0.387[Gr^Pr/[l+( .559/Pr) 9/16 ] 16/9 ^ 1/6 j 2 
where: ^ , Pr, and Gr are based on the fluid 

temperature , 

Tf = (T5 + 
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and Gr^ is based on the maximum diameter of the airship (see 
the temperature profile, Figure 4) . 

Under forced convection conditions, the airship 
external surface is idealized as a simple flat plate with 
the equivalent length of the airship and the area of the 
outer envelope. Equations 5-79 and 5-83 of reference 8 are 
manipulated to yield the average equivalent convective heat 
transfer coefficient: 

h = (k/L)ReLPr^/^{ [ .2275/ (log ReL)^-^®'^] 

- [850/ReL]} 
assumptions are: 

- laminar-turbulent boundary layer transition occurs at 
Re= 5x10^ 

- boundary layer is attached over the entire surface 

- the additional three dimensional body of rotation 
effects at the nose and tail are neglected 

- isothermal surface 

The derived convective heat transfer coefficients 
are then substituted into the linear heat transfer equation: 
q = hA(T 2 - . 

3 . Free Convection Between the Envelopes 

Here, due to the limitations of empirical data, the 
rough idealization of concentric horizontal cylinders is 
used. Equations 7-60 through 7-63, again from reference 8, 
are the closest approximation to the physical model 
available. These equations are based on isothermal 
surfaces, and a constant spacing in the annulus between the 
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Figure 4: Temperature Profile 



cylinders. A major problem here is that this empirical data 
is based on GrPr factors up to only 10®. The airship model 
GrPr combination is on the order of 10^^. However, since 
GrPr = 10° IS well into a turbulent convective boundary 
layer, it is assumed there are no significant changes up to 
GrPr = 10^^. Then: 

kg = k(0.4) [Gr^Pr]°-2 

where: Gr^ is based on the spacing between the cylinders, 
^out - ^in = ^ 

and the properties ^ , k are evaluated at T^ 

Tj^ = the "volume mean" temperature 

Rm = (Ri„ + Rout>/2 
then: q = 2 //Lkg(T 3 - T^)/ In 

4 . Free Convection Inside the Inner Envelope 

Modelling the convection in the inner envelope 
proves to be the most elusive problem. Empirical equations 
based on heat transfer from a uniform temperature fluid to a 
container wall are not adequate due to the additional 
convective turbulence in this volume gained from whatever 
heat source is used. Idealizing a heating duct and the 
envelope as concentric bodies results in Grashof numbers 
several orders of magnitude higher even than that between 
the envelopes. It is believed that completely different 
convective patterns will result, and that the application of 



28 



concentric body equations to this physical model is 
questionably accurate. 

The most applicable heat transfer coefficient 
available for this model can be derived from Stefan's 
article on hot air balloons (reference 5) . Although the 
convective activity for the airship heater or heat source is 
less violent than the open flame of the hot air balloon, the 
scale of convective activity is of the same order, and the 
heat transfer mechanisms at the envelope wall are assumed to 
be similar. Accordingly, Stefan's convective heat transfer 
coefficient, based on hot air balloon flight data, is 
adapted for the Helium in the airship by assuming 
Nu = (constant) (GrPr)^ 

as in the standard empirical forms. The exponent, n, is 
approximated as 0.25, similar to the value used in the 
volume between the envelopes. Then the constant is 
cancelled by: 

N%e = NUaij- (GrPr)H3°*25 / (GrPr) ^ * 25 

and 

^He = ^He ^air [ (^He/-^He f air) 2 ] 0 • 25 

From reference 4, 

hair [0-51 BTU/hr ff2 R] [Tq - T 2 ]^/^ 

and 

q = haij^(To - T 2 ) 

This approach results in values of h double that 
derived from the application of the standard empirical 
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equations previously referenced. It is believed to be most 
accurate due to the similarity of convective heating action 
and size scale. 

Additional complications arise in this internal 
volume if a steam injection method is used as a heating 
source, as described in section II, B of this paper. Basic 
heat transfer texts such as references 8 and 10 discuss this 
difficult problem at length. Convective heat transfer will 
decrease due to the reduced temperature of wet steam from 
that of a simple liquid fuel heater, and a small insulative 
effect of the condensate film on the envelope is gained. 

The significant change, however, is seen when the heat of 
vaporization is transferred almost directly to the envelope 
wall during condensation. A direct estimate of the heat 
transfer characteristics is not attempted for this case, but 
data based on opposite extremes is obtained in section IV. 
The optimistic extreme is assumed as no increase in the heat 
transfer coefficient. The pessimistic extreme is the case 
where the envelope wall temperature is the same as mean 
fluid temperature (or an infinite heat transfer 
coefficient) . 

C. RADIATION HEAT TRANSFER 
1 . General 

Initially, the low temperature expected in this 
airship would lead one to neglect the radiation heat 
transfer. The insulative effect of the inter-envelope air 
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space, however, reduces the convective heat loss to the 
order of the radiative heat loss between the envelopes for 
expected envelope emissivity values. 

Envelope emissivity, then, becomes of primary 
importance. Developments in aerospace metallic coatings 
used to control temperature of satellites provides the 
promise of desired characteristics. Reference 11 discusses 
vacuum metallizing or solution processing techniques of 
depositing a thin metal coating onto mylar film which will 
result in an opaque metal thickness of 5 x 10~^ inches. 

This thickness results in about 25 pounds of Aluminum for 
each envelope surface coated, and emissivities less than 
0 . 1 . 

The addition of the non-linear radiative heat 
transfer term in a set of linear equations becomes a major 
concern. Section D describes how the problem is handled in 
this model, and the problem of iteration convergence. 

2 . Radiation to the Inner Envelope 

The straightforward solution of the radiative heat 
transfer between an enclosed body and the enclosing surface 
leads to equation 8-43 of reference 8: 

<3r = O' ^htr(Thtr'*-'^2'‘) 

<Vfhtr> + [*htr/''llIV£i - H 

where and^j^^^ can be chosen appropriately. In 

the results section, from program HOTSHIP, these values are 
used pessimistically as: 



31 



= 0.1 (typical of metal surface) 
and results are derived for a range of envelope 
emissivities . 

Computations are more complex for steam injection. 
Strong absorption lines for the water vapor lead to 
calculations of radiative heat transfer based on the mean 
free path computations as shown in chapter 8 of reference 8, 
and in reference 6. Due to the absence of the heater with a 
high surface temperature, however, the total radiative loss 
is expected to be less with the steam injection. This 
computation is not completed in this model, and program 
HOTSHIP assumes the radiative heat loss to be the same for 
both heating methods. 

3 . Radiation Between the Envelopes 

The greatest radiative heat loss problem is between 
the envelopes, due to the large radiative surface areas. 
Again, the enclosed body equation is used: 

= <r A^n(T3'»-T4'‘) 

(l/£in) + [Ain/*oHV£o ‘ H 

and results are obtained for envelope emissivities from 0 to 

0.9. 

4 . Radiative Loss from the Outer Envelope 
References 5 and 6 have a detailed treatment of 

radiation heat transfer between a hot balloon surface and 



the surrounding environment. A direct application to this 
airship model can be made, but this is not completed for 
three reasons. First, from the high convective heat loss 
only, the external envelope temperature is normally within 
two degrees F of the outside air temperature. Second, the 
expected operational conditions of this airship are assumed 
to vary to extremes, and results from reference 5 show that 
radiative heat loss is highly variable in these conditions. 
Accounting for these variations is not considered 
appropriate or necessary in this initial model. Third, due 
to expected night and foul weather operations, the external 
envelope emissivity should be small, minimizing the external 
radiative heat loss. Under specified environmental 
conditions, program HOTSHIP can be modified to include this 
external radiation term. 

D. THE RESULTING SET OF LINEARIZED EQUATIONS 

Discretizing the temperature profile as in Figure 4, 
then, the linearized set of heat transfer equations becomes 
a linear system of five equations in five unknowns for this 
particular physical model. 

At any surface, 

g = ‘^conv + ^r 
In the inner envelope, then, 

g = Iconv + <ar 

= h^AiCTo - Tj) + O- Ahtr<'^htr''-'^2'') 

(V£htr) + [Ahtr/'^iHV£i - 1] 
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and to linearize, an initial value is substituted for T2 in 
the qj, term, and q^. is treated as a constant in the linear 
system. Then, 

q = hj^Aj^[TQ - T2) + QrO-2 
Similarly, between the envelopes, 

q = 2 a L(T3 - T^) + <r Ahtr(Thtr"^-T2^) 

Iri (^out/^in^ htr^ ^^htr^^i^ ^ i ~ 

= 2 /7 kg L(T3 - T4) + ' 3 r 3-4 

ln(Rout/I^in) 

where initial values of T3 and are substituted to make 
qr3_4 a constant in the system. 

Rearranging terms, the resulting set of equations 
becomes : 



T 



2 



+ 



q/hiAi = Tq + (q3.o_2)/hj^Ai 



T 



2 



-q = 0 



T-, - T. _ - g In (Rout/Rin) =-g„o . In (Rout/Rin) 

2 ;7 kgL 2)7 kgL 

T4 - T5 - q to/ko^o = 0 

'^5 “ q/^ 4 Ao = '^inf 

In matrix form. 



"1000 1/h^A^ 




r s 

^2 


1 


To + (qr0-2^/^l^i 


1-10 0 -tj^/k^Aj^ 




T 3 




0 


01-10 -InfRout/Rin) 




Td 


= 


-q^-, . In fRout/Rin) 


2 77 kgL 




4 




2 7! kgL 


001-1 -tg/kgAg 




T 5 




0 


0001 -l/h^Ag 

y 




q 




■^inf 
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SECTION 5: VALUES OBTAINED FOR INCREASE ENVELOPE HEAT LOSS 

RESISTANCE 

FOR A FACTOR OF 2 INCREASE: 

1 

INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 

.1, .9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 

■> 

50,540 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 8 
DQDTN = 39.79 DMDT = 6.99 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.58 PDLTOT= 3.98 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

550.00 548.83 548.47 540.41 540.16 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.032281 0.008935 0.553427 0.006463 0.003923 

DQDTRO = 3.54 DQDTR3 = 25.23 DQDTN = 39.79 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 9 
DQDTN = 86.46 DMDT = 15.18 

RIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

PDL = 5.14 PDLTOT= 9.51- 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

560.00 557.82 557.06 540.90 540.34 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.026212 0.008855 0.474857 0.006463 0.003923 

DQDTRO = 3.40 DQDTR3 = 52.43 DQDTN = 86.46 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR 0 9 



DQDTN = 137.58 DMDT = 

RIN = 46.68 AIN = 

VIN = 2192307.69 DVH20 = 

PDL = 7.68 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

570.00 566.89 565.68 

THE MATRIX VALUES A,B,C,D,E 
0.023150 0.008776 0.431857 

DQDTRO = 3.24 DQDTR3 = 



24.16 

104666.37 

121494.44 

14.32 

,T2,T3,T4,T5) IS: 

541.43 540.54 

ARE: 

0.006463 0.003923 

81.42 DQDTN = 137.58 



NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 10 
DQDTN = 192.52 DMDT = 33.81 

RIN = 47.15 AIN = 105586.99 

VIN = 2230769.23 DVH20 = 192264.72 

PDL = 10.19 PDLTOT= 20.64 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

580.00 576.00 574.32 542.00 540.76 
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THE MATRIX VALUES A,B,C,D,E ARE: 

0.020918 0.004350 0.399194 0.003231 0.003923 

DQDTRO = 3.08 DQDTR3 = 116.51 DQDTN = 200.58 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 
DQDTN = 262.07 DMDT = 46.02 

RIN = 47.61 AIN = 106500.01 
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VIN ~= 2269230.77 DVH20 = 333175.90 

PDL = 12.65 PDLTOT= 30.69 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS; 

590.00 584.95 583.82 541.87 541.03 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.019469 0.004312 0.376496 0.003231 0.003923 

DQDTRO = 2.92 DQDTR3 = 150.64 DQDTN = 262.07 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 327.36 DMDT = 57.49 

RIN = 48.07 AIN = 107405.63 

VIN = 2307692.31 DVH20 = 453004.31 

PDL = 15.10 PDLTOT= 39.47 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

6^.00 594.05 592.65 542.34 541.28 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.018332 0.004276 0.358000 0.003231 0.003923 

DQDTRO = 2.74 DQDTR3 = 186.84 DQDTN = 327.36 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 396.42 DMDT = 69.61 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.51 PDLTOT= 52.59 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

610.00 603.15 601.47 542.84 541.56 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.017400 0.004241 0.342319 0.003231 0.003923 

DQDTRO = 2.55 DQDTR3 = 225.15 DQDTN = 396.42 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 469.23 DMDT = 82.40 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.89 PDLTOT= 76.87 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

620.00 612.24 610.27 543.36 541.84 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016613 0.004206 0.328650 0.003231 0.003923 

DQDTRO = 2.36 DQDTR3 = 265.63 DQDTN = 469.23 
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VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 21.51 PDLTOT= 100.51 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

630.00 589.08 571.00 571.00 557.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.009450 0.004172 0.000000 0.003231 0.003923 

DQDTRO = 2.84 DQDTR3 = 0.00 DQDTN = 4333.20 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 9 
DQDTN = 4897.97 DMDT = 860.13 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVK20 = 2272776.29 

PDL = 23.82 PDLTOT= 144.22 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

640.00 595.32 575.04 575.04 559.21 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.009128 0.004139 0.000000 0.003231 0.003923 

DQDTRO = 2.71 DQDTR3 = 0.00 DQDTN = 4897.97 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 9 
DQDTN = 5472.35 DMDT = 961.00 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.10 PDLTOT= 258.81 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

650.00 601.63 579.15 579.15 561.47 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.008844 0.004107 0.000000 0.003231 _0. 003923 

DQDTRO = 2.59 DQDTR3 = 0.00 DQDTN = 5472.35 

NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 9 
DQDTN = 6055.73 DMDT = 1063.44 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 28.36 PDLTOT= 452.45 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

660.00 608.00 583.32 583.32 563.76 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.008590 0.004075 0.000000 '0.003231 0.003923 

DQDTRO = 2.45 DQDTR3 = 0.00 DQDTN = 6055.73 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 9 
DQDTN = 6647.57 DMDT = 1167.38 
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SECTION 3s VALUES OBTAINED FOR T3 = T4 



INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 

7 

.1, .9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 

7 

50,540 

NUMBER OF ITERATIONS AREsRIN 5 ,FOR Q 7 
DQDTN = 332.86 DMDT = 58.45 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.32 PDLTOT= 3.74 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) ISs 

550.00 543.87 542.38 542.38 541.31 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.018624 0.004468 0.000000 0.003231 0.003923 

DQDTRO = 3.62 DQDTR3 = 0.00 DQDTN = 332.86 

NUMBER OF ITERATIONS AREsRIN 5 ,FOR Q 7 
DQDTN = 751.63 DMDT = 131.99 

RIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

PDL = 4.73 PDLTOT= 9.13 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) ISs 

560.00 548.71 545.38 545.38 542.95 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.015099 0.004427 0.000000 0.003231 0.003923 

DQDTRO = 3.54 DQDTR3 = 0.00 DQDTN = 751.63 

NUMBER OF ITERATIONS ARE ; RIN 5 , FOR Q 8 
DQDTN = 1207.26 DMDT = 212.01 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.16 PDLTOT= 13.86 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

570.00 553.93 548.64 ' 548.64 544.74 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.013346 0.004388 0.000000 0.003231 0.003923 

DQDTRO = 3.46 DQDTR3 = 0.00 DQDTN = 1207.26 



NUMBER OF ITERATIONS AREsRIN 5 ,FOR Q 8 
DQDTN = 1688.06 DMDT = 296.44 

RIN = 47.15 AIN = 105586.99 

VIN = 2230769.23 DVH20 = 192264.72 

PDL = 9.60 PDLTOT= 20.14 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) ISs 

580.00 559.42 552.08 552.08 546.62 
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RIN 




51.22 


AIN = 


113552.84 




VIN 


= 


2576923.08 


DVH20 = 


35853817.74 




PDL 




31.99 


PDLTOT= 


1859.37 




THE 


TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) 


IS; 


670 


.00 


670.00 


665.92 


547.22 


543.96 



THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004045 0.272200 0.003231 0.003923 

DQDTRO = 0.92 DQDTR3 = 573.25 DQDTN = 1009.32 

NUMBER OF ITERATIONS AREtRIN 5 ,FOR Q 6 
DQDTN = 1119.72 DMDT = 196.63 

RIN = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVK20 = **-/^*^^-^*^** 

PDL = 34.20 PDLTOT= -799.89 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS; 

680.00 630.00 675.50 548.01 544.39 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004014 0.264220 0.003231 0.003923 

DQDTRO = 0.62 DQDTR3 = 637.19 DQDTN = 1119.72 
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NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 6 
DQDTN = 616.34 DMDT = 108.24 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.79 PDLTOT= 100.52 

THE TEMPERATURE PROFILE (T0,T2,T3 ,T4,T5) IS: 

630.00 630.00 627.43 544.41 542.42 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004172 0.310342 0.003231 0.003923 

DQDTRO = 1.96 DQDTR3 = 348.83 DQDTN = 616.34 

NUMBER OF ITERATIONS AREtRIN 4 ,FOR Q 6 
DQDTN = 707.53 DMDT = 124.25 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 25.15 PDLTOT= 143.52 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

640.00 640.00 637.07 545.06 542.78 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004139 0.299614 0.003231 0.003923 

DQDTRO = 1.71 DQDTR3 = 400.44 DQDTN = 707.53 

NUMBER OF ITERATIONS ARE:RIN 4 .FOR Q 6 
DQDTN = 803.35 DMDT = 141.08 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 27.46 PDLTOT= 256.09 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

650.00 650.00 646.70 545.75 543.15 

THE MATRIX VALUES A,B,C.D,E ARE: 

0.000000 0.004107 0.289787 0.003231 0.003923 

DQDTRO = 1.46 DQDTR3 = 454.98 DQDTN = 803.35 

NUMBER OF ITERATIONS ARE.-RIN 5 .FOR Q 6 
DQDTN = 903.90 DMDT = 158.73 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.74 PDLT0T= • 446.08 

THE TEMPERATURE PROFILE (T0.T2.T3.T4.T5) IS: 

660.00 660.00 656.32 546.47 543.55 

THE MATRIX VALUES A.B.C.D.E ARE: 

0.000000 0.004075 0.280690 0.003231 0.003923 

DQDTRO = 1.19 DQDTR3 = 512.55 DQDTN = 903.90 

NUMBER OF ITERATIONS ARE: RIN 5 .FOR Q 6 
DQDTN = 1009.32 DMDT = 177.25 
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SECTION 2: VALUES OBTAINED FOR T2 = TO 



INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 

7 

.1,.9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 

7 

50,540 



NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 5 
DQDTN = 47.35 DMDT = 8.32 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.67 PDLTOT= 4.07 

THE TEMPERATURE PROFILE (T0,T2 ,T3 ,T4 ,T5) IS; 

550.00 550.00 549.79 540.34 540.19 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.000000 0.004468 0.535970 0.003231 0.003923 

DQDTRO = 3.52 DQDTR3 = 29.72 DQDTN = 47.35 

NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 5 
DQDTN = 102.26 DMDT = 17.96 

RIN = 46.22 AIN = 103737.93 

VIN = 2153846.15 DVH20 = 79283.17 

PDL = 5.31 PDLTOT= 9.67 

THE TEMPERATURE PROFILE (TO ,T2,T3 ,T4,T5) IS; 

560.00 560.00 559.55 540.73 540.40 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004427 0.460701 0.003231 0.003923 

DQDTRO = 3.36 DQDTR3 = 61.42 DQDTN = 102.26 

NUMBER OF ITERATIONS ARE:RIN 5 ,F0R Q 5 
DQDTN = 162.32 DMDT = 28.51 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.92 PDLTOT= 14.54 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

570.00 570.00 569.29 ' 541.16 540.64 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004388 0.419335 0.003231 0.003923 

DQDTRO = 3.19 DQDTR3 = 95.25 DQDTN = 162.32 



NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 6 
DQDTN = 227.01 DMDT = 39.87 

RIN = 47.15 AIN = 105586.99 

VIN = 2230769.23 DVH20 = 192264.72 

PDL = 10.49 PDLTOT= 20.90 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

580.00 580.00 579.01 541.62 540.89 
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THE MATRIX VALUES A,B,C,D,E ARE: 

0.013889 0.004045 0.277434 0.003231 0.003923 

DQDTRO = 1.26 DQDTR3 = 503.05 DQDTN = 891.11 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 12 
DQDTN = 987.43 DMDT = 173.40 

RIN = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVH20 = *********** 

PDL = 33.48 PDLTOT= -809.11 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

680.00 666.70 662.73 547.06 543.87 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.013485 0.004014 0.269284 0.003231 0.003923 

DQDTRO = 1.01 DQDTR3 = 557.89 DQDTN = 987.43 
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vapor at saturation level in the Helium becomes more and 
more rapid as temperature is increased. As previously 
noted, the loss of insulative space from the water vapor 
volume is not accounted for in program HOTSHIP, and the heat 
loss rates shown are increasingly optimistic with 
temperature rise. This effect is not more than a few 
percent at 120 deg F (delta T = 60) but by 160 deg F, the 
inner volume would be greater than the outer volume, 
implying that by this point the inter envelope insulative 
effect is absent, and the heat flux will actually be many 
times higher as shown previously in this section. 

If losses due to condensation on the inner envelope are 
much higher than expected, the pessimistic approximation of 
T 2 = Tq is made by setting the matrix value of "A" to zero. 
Results in Appendix B show only a 12.9% increase in heat 
loss for this condition. 

B. ESTIMATES OF ACCURACY 

The best accuracy is obtained from the conduction 
equations in the envelope materials. For the material 
described, the conductivity is estimated to be within 20% of 
the average of that to be obtained in an actual envelope. 

The one dimensional conduction equation is completely 
applicable. Unfortunately, this is among the least critical 
conditions due to the ineffective insulation of the 
envelopes. Even if the conductivity of the envelopes was 
halved, the effect on heat loss is minor. In Appendix B, 
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Figure 9: Change in Lift vs Q 



Although some heat loss resistance is exhibited at the 
boundary layers on the inside and outside, the results file 
in Appendix B clearly shows that the inter-envelope volume 
provides most of the insulative effects. More than 8/10 of 
the temperature drop is seen across this air gap under 
"normal" cruise conditions, implying that the heat loss 
would be several times higher in a conventional airship 
design. In fact, if the matrix value "C" is set to zero in 
program HOTSHIP, the heat flux is seen to increase by a 
factor of 8.4 (see Appendix B) . This result does not 
account for the additional radiative heat loss which will 
occur due to the significant temperature increase of the 
outer envelope. 

This insulative gap is seen to be particularly critical 
as airship speed is increased from zero. Figure 10 shows 
that the heat flux quickly rises to an asymptotic value 
about twice that for zero velocity even in the concept 
design. This rise would be accentuated by the appropriate 
factor of more than eight without the inter envelope space. 
This implies that a standard airship moving at 50 knots will 
have an order of magnitude higher heat loss than a hot air 
balloon with the same shape and size, and 17 times the heat 
loss of the concept design at zero velocity. 

Figures 8 and 9 show a dramatic difference in lift 
increase (for the same heat input) when the steam injection 
heating concept is used. The volume increase due to water 
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C FORCED CONVECTION PARAMETERS OUTSIDE THE AIRSHIP, BASED ON A 
C CRITICAL REYNOLDS NUMBER RECRIT=500000 
RHOI=0 . 0022846*540 .0/TI 

MUI=((4.3339E-7-3.85603E-7)/90.0)*(TI-540.0)+3.85605E-7 
REL=RHOI*UI*6076 . 0*L/ (MUI*3600 .0) 
PRI=((.697-.708)/90.0)*(TI-540.0)+.708 

KINF=( ( (.0173513-. 0151615)/90.0)*(TI-540.0)+.0151615)/3600.0 

H4= ( (0 . 2275/ ( (LOGIO (REL) ) **2 . 584) ) -850 . 0/REL) *KINF*REL* 
+(PRI**0.333333333)/L 

C SETTING THE INITIAL VALUES FOR THE LINEAR SYSTEM ITE.RATION 

400 T2=T0-(T0-TI)/5.0 

T3=T0-(T0-TI)/4.5 
T4=T0-(T0-TI)*4. 0/5.0 
T5=T0-(T0-TI)*4. 5/5.0 
DQDT=480 
EHTR=0 . 1 
AHTR=500 
THTR=700 

J=1 

IF (UI.NE.0.0) GO TO 600 

C IF EXTERNAL FLOW IS FREE CONVECTION, COMPUTE THE HEAT TRANSFER 

C COEFFSCIENT FOR UI=0, IF NOT, H4 IS ALREADY COMPUTED 

500 TFL=(T5+TI)/2 

PRFL=((.697-.708)/90.0)*(TFL-540.0)+.708 

KFL=(( (.0173513-. 0151615)/90. 0)*(TFL-540. 0)+. 0151615 )/3600.0 
MUFL=((4.3339E-7-3.85605E-7)/90.0)*(TFL-540.0)i-3.85605E-7 
RH0FL=0 . 0022846*540 . 0/TFL 

GRFL=32 . 174* (T5-TI ) *8 . 0* (ROUT**3 ) * (RH0FL**2 ) / (TFL* (MUFL**2) ) 

H4=(KFL/L)*( (.6+. 387*(GRFL*PRFL/(1.0+(.5’59/PRFL)**. 5625) **1.7778) 
+** .1666667 )**2) 

C FREE CONVECTION BETWEEN THE ENVELOPES BASED ON HORIZONTAL 
C CYLINDRICAL ANNULUS, ISOTHE.R.MAL ENVELOPES, AND CONSTANT 

C ENVELOPE SPACING 

600 RMEAN=(ROUT+RIN)/2.0 

TMEAN= ( ( (RMEAN**3 ) - (RIN**3 ) ) *T3+ ( ( ROUT**3 ) - ( RMEAN**3 ) ) *T4 ) / 

+ ( ( ROUT**3 ) - ( RIN**3 ) ) 

RHOM=0 .0022846*540 . 0/TMEAN 

MUMN=3.85605E-7+((4.3339E-7-3.85605E-7)/90.0)*(TMEAN-540.0) 

PRMN=0.708+((.697-.708)/90.0)*(TMEAN-540.0) 

KMEAN=(0.0151615+( (.0173513-. 0151615) /90.0)*(TME AN-540. 0))/3600.0 
GRD=32.174*(T3-T4)*((ROUT-RIN)**3)*(RHOMN**2)/(TMEAN*(MUMN**2)) 
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THE MATRIX VALUES A,B,C,D,E ARE: 

0.021137 0.008700 0.402213 0.006463 0.003923 

DQDTRO = 3.08 DQDTR3 = 112.16 DQDTN = 192.52 



NUMBER OF ITERATIONS ARE:RIN 5 ,F0R Q 10 
DQDTN = 251.11 DMDT = 44.10 

RIN = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 12.66 PDLTOT= 30.70 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5 ) IS: 

590.00 585.12 582.95 542.61 540.99 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.019680 0.008625 0.379454 0.006463 0.003923 

DQDTRO = 2.91 DQDTR3 = 144.80 DQDTN = 251.11 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 313.19 DMDT = 55.00 

RIN = 48.07 AIN = 107405.63 

VIN = 2307692.31 DVH20 = 453004.31 

PDL = 15.11 PDLTOT= 39.48 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

600.00 594724 591.57 543.15 541.23 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.018538 0.008552 0.360910 0.006463 0.003923 

DQDTRO = 2.74 DQDTR3 = 179.32 DQDTN = 313.19 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 378.69 DMDT = 66.50 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.52 PDLTOT= 52.59 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS : 

610.00 603.38 600.17 543.93 541.49 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.017602 0.008481 0.345190 0.006463 0.003923 

DQDTRO = 2.55 DQDTR3 = 215.78 DQDTN = 378.69 

NUMBER OF ITERATIONS ARE:RIN. 5 ,FOR Q 11 
DQDTN = 447.65 DMDT = 78.61 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.91 PDLTOT= 76.87 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

620.00 612.51 608.75 544.65 541.76 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.016819 0.008412 0.331477 0.006463 0.003923 

DQDTRO = 2.36 DQDTR3 = 254.29 DQDTN = 447.65 



NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 
DQDTN = 520.00 DMDT = 91.32 

RIN = 49.44 AIN = 110079.94 
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VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.25 PDLTOT= 100.51 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

630.00 621.64 617.31 545.40 542.04 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016136 0.008344 0.319288 0.006463 0.003923 
DQDTRO = 2.15 DQDTR3 = 294.80 DQDTN = 520.00 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 
DQDTN = 595.80 DMDT = 104.63 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.57 PDLTOT= 143.82 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

640.00 630.77 625.84 546.19 542.34 

T.HE MATRIX VALUES A,B,C,D,E ARE: 

0.015535 0.008278 0.308269 0.006463 0.003923 

DQDTRO = 1.94 DQDTR3= 337.41 DQDTN = 595.80 

NUMBER OF ITERATIONS ARE;RIN 4 ,FOR Q 11 
DQDTN = 675.08 DMDT = 118.55 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVK20 = 4419910.85 

PDL = 26.85 PDLTOT= 257.31 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

650.00 639.90 634.35 547.01 542.65 

JHE MATRIX VALUES A,B,C,D,E ARE: 

0.015000 0.008214 0.298185 0.006463 0.003923 

DQDTRO = 1.72 DQDTR3 = 382.16 DQDTN = 675.08 

NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 11 
DQDTN = 757.89 DMDT = 133.09 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.10 PDLTOT= 449.04 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

660.00 649.02 642.84 547.87 542.97 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014518 0.008151 0.288857 0.006463 0.003923 

DQDTRO = 1.49 DQDTR3 = 429.12 DQDTN = 757.89 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 12 

DQDTN = 844.26 DMDT = 148.26 

RIN = 51.22 AIN = 113552.84 

VIN = 2576923.08 DVH20 = 35853817.74 

PDL = 31.32 PDLTOT= 1875.65 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

670.00 658.13 651.30 548.77 543.31 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014076 0.008089 0.280164 0.006463 0.003923 

DQDTRO = 1.24 DQDTR3 = 478.27 DQDTN = 844.26 
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NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 12 
DQDTN = 934.32 DMDT = 164.08 

RIN = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVH20 = *********** 

PDL = 33.50 PDLTOT= -808.84 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

680.00 667.24 659.74 549.70 543.67 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.013674 0.008029 0.271992 0.006463 0.003923 

DQDTRO = 0.99 DQDTR3 = 529.77 DQDTN = 934.32 



FOR A FACTOR OF 10 INCREASE: 



INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 
•> 

.1, .9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 

•> 

50,540 

NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 9 
DQDTN = 31.55 DMDT = 5.54 

RIN = 45.75 AIN = 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.60 PDLTOT= 4.00 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

550.00 549.03 547.62 541.14 540.12 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.034562 0.044676 0.577787 0.032314 0.003923 

DQDTRO = 3.54 DQDTR3 = 20.33 DQDTN = 31.55 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR 0 9 



DQDTN = 67.25 DMDT = 

RIN = 46.22 AIN = 

VIN = 2153846.15 DVK20 = 

PDL = 5.17 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

560.00 558.21 555.23 

THE MATRIX VALUES A,B,C,D,E 
0.028005 0.044273 0.497466 

DQDTRO = 3.39 DQDTR3 = 



11.81 

103737.93 

79283.17 

9.54 

,T2,T3,T4,T5) IS: 

542.44 540.26 

ARE: 

0.032314 0.003923 

41.52 DQDTN = 67.25 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 105.57 DMDT = 18.54 

RIN = 46.68 AIN = 104666.37 



VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.72 PDLTOT= 14.36 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

570.00 567.47 562.83 543.83 540.41 
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NUMBER OF ITERATIONS AREiRIN 5 ,FOR Q 12 
DQDTN = 1013.73 DMDT = 178.02 

RIN = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVH20 = *********** 

PDL = 33.21 PDLTOT= -812.60 

THE TEMPERATURE PROFILE (TO ,T2,T3 ,T4 ,T5) IS; 

680.00 666.43 662.36 543.28 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.013397 0.004014 0.267670 0.003231 0.000000 

DQDTRO = 1.02 DQDTR3 = 568.82 DQDTN = 1013.73 
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THE MATRIX VALUES A,B,C,D,E ARE: 

0.024750 0.043880 0.453521 0.032314 0.003923 

DQDTRO = 3.23 DQDTR3 = 63.66 DQDTN = 105.57 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR 0 10 



DQDTN = 146.21 DMDT = 

RIN = 47.15 AIN = 

VIN = 2230769.23 DVH20 = 

PDL = 10.24 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

580.00 576.76 570.40 

THE MATRIX VALUES A,B,C,D,E 
0.022670 0.043498 0.423124 

DQDTRO = 3.07 D0DTR3 = 



25.68 

105586.99 

192264.72 

20.69 

,T2,T3,T4,T5) IS: 
545.30 540.57 



ARE: 

0.032314 0.003923 

86.89 DQDTN = 146.21 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR 0 11 



DQDTN = 188.96 DMDT = 

RIN = 47.61 AIN = 

VIN = 2269230.77 DVH20 = 

PDL = 12.73 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

590.00 586.06 577.91 

THE MATRIX VALUES A,B,C,D,E 
0.021164 0.043125 0.399790 

DQDTRO = 2.90 DODTR3 = 



33.18 

106500.01 

333175.90 

30.75 

,T2,T3,T4,T5) IS: 

546.85 540.74 

ARE; 

0.032314 0.003923 

111.26 DQDTN = 188.96 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR 0 11 



DQDTN = 233.65 DMDT = 

RIN = 48.07 AIN = 

VIN = 2307692.31 DVH20 = 

PDL = 15.18 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

600.00 595.39 585.40 

THE MATRIX VALUES A,B,C,D,E 
0.019974 0.042761 0.380820 

DQDTRO = 2.71 DQDTR3 =. 



41.03 

107405.63 

453004.31 

39.53 

,T2,T3,T4,T5) IS: 

548.47 540.92 

ARE: 

0.032314 0.003923 

136.68 DQDTN = 233.65 



NUMBER OF ITERATIONS ARE;RIN. 5 ,FOR Q 11 
DQDTN = 280.24 DMDT = 49.21 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.60 PDLTOT= 52.64 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5 ) IS: 

610.00 604.72 592.84 550.16 541.10 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.019002 0.042406 0.364745 0.032314 0.003923 

DQDTRO = 2.52 DQDTR3 = 163.22 DQDTN = 280.24 



NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 11 
DQDTN = 328.70 DMDT = 57.72 

RIN = 48.98 AIN = 109195.42 
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VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.99 PDLTOT= 76.90 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

620.00 614.07 600.24 551.91 541.29 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.018182 0.042060 0.350738 0.032314 0.003923 

DQDTRO = 2.32 DQDTR3 = 190.90 DQDTN = 328.70 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 12 
DQDTN = 378.94 DMDT = 66.55 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.35 PDLTOT= 100.51 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

630.00 623.42 607.61 553.73 541.49 

THE MATRIX VALUES A,B,C,D,E ARE: 

0 .017 4 70 0.041 722 _0 . 3 3 82 9 6 0.032314 0.003923 

DQDTRO = 2.11 DQDTR3 = 219.69 DQDTN = 378.94 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 12 
DQDTN = 431.05 DMDT = 75.70 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.67 PDLTOT= 143.77 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

640.00 632.77 614.93 555.62 541.69 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.016849 0.041392 0.327034 0.032314 0.003923 

DQDTRO = 1.89 DQDTR3 = 249.70 DQDTN = 431.05 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 12 
DQDTN = 484.96 DMDT = 85.16 

RIN = ' ~50.33 AIN = 111828.'^ 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.95 PDLTOT= 257.11 

THE TEMPERATURE PROFILE (TO ,T2,T3 ,T4,T5) IS: 

650.00 642.12 622.21 557.57 541.90 

THE MATRIX VALUES A,B,C,D,E ARE-* 

0.016297 0.041070 0.316728 0.032314 0.003923 

DQDTRO = 1.66 DQDTR3 = 280.90 DQDTN = 484.96 

NUMBER OF ITERATIONS ARE-.RIN 5 ,FOR Q 12 
DQDTN = 540.69 DMDT = 94.95 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.21 PDLTOT= 448.55 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

660.00 651.48 629.44 559.59 542.12 

THE MATRIX VALUES A,B,C,D,E ARE : 

0.015799 0.040754 0.307195 0.032314 0.003923 

DQDTRO = 1.42 DQDTR3 = 313.30 DQDTN = 540.69 
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APPENDIX A 



PROGRAM HOTSHIP 

The following FORTRAN program, as described in the text, 
is used to iterate the linearized heat transfer system of 
equations. Program inputs are air temperature, and 

envelope emissivities . Several outputs, including the heat 
flux, q, and temperature profile values, are written to the 
screen and to files for use with the graphics program used 
to compile figures 5 through 10. 
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C THIS PROGRAM ITERATES A LINEAR HEAT TRANSFER MODEL OF AN AIRSHIP 

C TO SOLVE FOR THE TEMPERATURE PROFILE AND THE RATE OF HEAT TRANSFER 

REAL MUI , MUMN , HI , H4 , KMEAN , KINF , KEQV , KIN , KOUT , GRD , KFL , MUFL , GRFL 
REAL KHE , KAI R , MUAI R , MUHE , LH20 , MH20 , HI AI R 

C INPUTS ARE; THE ENVELOPE EMISSIVITIES , VELOCITY, AND AIR TEMP 
PI=3. 1415927 

10 PRINT*, 'DO YOU WANT TO INCLUDE THE RADIATION HEAT TRANSFER' 

PRINT* , ' EFFECTS ? ( 1=YES , 2=NO ) ' 

READ*, Y 

IF (Y.EQ.2.0) GO TO 20 

PRINT*, 'INPUT THE EMISSIVITIES OF THE INNER AND OUTER' 

PRINT*, 'ENVELOPES, EIN AND EOUT' 

READ*, EIN, EOUT 

20 PRINT *, 'INPUT U INFINITY (UI,KTS), AND T INF (DEG R) ' 

READ *, UI,TI 

C COMPUTING THE DIMENSIONS OF THE INNER ENVELOPE 
C WITH THESE GIVEN OUTER ENVELOPE DIMENSIONS 
VTOT=4.0E6 
ROUT=66.5 
AOUT=1.4213E5 
L=398.9 

C THIS LOOP INCREMENTS THE HELIUM TEMP IN 10 DEG STEPS 
DO 1200 K=l,14 
T0=TI+10.0*K 
M=1 

VIN=2.0E6*T0/520 
RINIT=50 . 0 

100 RTRY=( -3 . 0*VIN/ ( 5 . 0*PI* (RINIT-3 . 0*L/5 . 0 ) ) ) **0 . 5 

IF (ABS(RTRY-RINIT).LE. 0.001) GO TO 200 
RINIT=RTRY 

IF (M.GE.20) GO TO 200 
M=M+1 
GO TO 100 
200 RIN=RTRY 

AIN= 2 . 0*PI*RIN**2+2 . 0*PI*RIN* (L-3 . 0*RIN)+PI* (RIN**2 ) *2 . 236068 

C SPECIFYING THE PARAMETERS OF THE ENVELOPE MATERIALS 

KIN=0. 03266/3600.0 
THIN=0. 05/12.0 
KOUT=KIN 
THOUT=THIN 

C IF UI=0, THEN FREE CONVECTION CONDITIONS EXIST 
IF (UI.EQ.0.0) GO TO 400 
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V. CONCLUSIONS 



For this relatively simple concept design, heat loss is 
small enough to allow large changes in lift by heating the 
lifting gas with engine waste or exhaust heat or small 
liquid fuel furnaces. For envelopes with low emissivities 
produced as discussed, these results, and Appendix C 
calculations show that more than a 30% increase might be 
expected by use of engine waste heat alone under moderate 
power conditions. If steam injection is used, a 100% 
increase might be expected, with a significant time and fuel 
penalty required for evaporating the required amount of 
steam. 

By using heated gas, vectored thrust, dynamic lift 
and/or other lift augmentation systems, the airship as 
conceived might achieve more than a 100% increase in lift, 
which will allow a tremendous operational flexibility in 
cargo load capacity, unrefueled range and endurance and 
mission equipment carried. Obvious advantages over rotary 
wing aircraft in load carrying capacity include a much 
larger weight capacity, low fuel useage, longer range, and 
no rotor blast effects on takeoff or landing. Advantages 
over fixed wing aircraft are vertical landing and takeoff 
and lower fuel useage. 
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IV. RESULTS 



A. PROGRAM OUTPUT DATA 

Figures 5 through 10 show the program HOTSHIP results 
graphically for the parameters of interest. All results 
shown are for ambient air temperature of 60 deg F. 

Certainly this airship concept will have the same response 
to atmospheric temperature and density changes as a standard 
airship. An additional minor effect will be the change in 
the inner envelope radius required to hold the same gross 
lift, which changes the insulative gap between the 
envelopes . 

It is quickly seen from Figures 5 through 7 that 
envelope emissivities will make a dramatic difference in 
heat loss. For organic materials, with emissivities greater 
than 0.9, the total heat loss rate is six or more times the 
convective heat loss alone ((£ = 0.0). Even with envelope 
emissivity of 0.1, the radiative and convective heat losses 
are on the same order. 

The curve for the specific results of 
£out “ shown, which is expected to correspond closely 

to the special case where the inner envelope only is vacuum 
metallized. Note that the heat loss for this case is less 
than if both envelopes have an emissivity of 0.2. 
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RIN = 51.22 AIN = 113552.84 

VIN = 2576923.08 DVH20 = 35853817.74 

PDL = 30.59 PDLT0T= 1893.36 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

670.00 614.45 587.56 587.56 566.08 

THE MATRIX VALUES A,B,C,D,E ARE ; 

0.008360 0.004045 0.000000 0.003231 0.003923 

DQDTRO = 2.31 DQDTR3 = 0.00 DQDTN = 6647.57 

NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 9 
DQDTN = 7247.43 DMDT = 1272.72 

RIN'. = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVH20 = *********** 

PDL = 32.78 PDLTOT= -818.05 

THE TEMPERATURE PROFILE (T0,T2.T3,T4,T5) IS; 

680.00 620.94 591.85 591.85 568.43 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.008151 0.004014 0.000000 0.003231 0.003923 

DQDTRO = 2.17 DQDTR3 = 0.00 DQDTN = 7247.43 
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THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

630.00 621.18 618.84 541.81 540.00 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.015844 0.004172 0.314878 0.003231 0.000000 

DQDTRO = 2.16 DQDTR3 = 314.35 DQDTN = 559.01 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 
DQDTN = 641.71 DMDT = 112.69 

RIN = 49.89 AIN = 110957.85 

VIN = 2461533.46 DVH20 = 2272776.29 

PDL = 24.35 PDLTOT= 143.94 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

640'; 00 630.25 627.59 542.07 540.00 

THE M.ATRIX VALUES A,B,C,D,E ARE: 

0.015248 0.004139 0.303883 0.003231 0.000000 

DQDTRO = 1.95 DQDTR3 = 360.30 DQDTN = 641.71 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 
DQDTN = 728.44 DMDT = 127.92 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.61 PDLTOT= 257.80 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

650.00 639.31 636.31 542.35 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014716 0.004107 0.293819 0.003231 0.000000 

DQDTRO = 1.73 DQDTR3 = 408.65 DQDTN = 728.44 

NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 11 
DQDTN = 819.31 DMDT = 143.88 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 28.84 PDLTOT= 450.24 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

660.00 648.36 645.02 542.65 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014236 0.004075 0.284508 0.003231 0.000000 

DQDTRO = 1.50 DQDTR3 = 459.49 DQDTN = 819.31 

NUMBER OF ITERATIONS ARE : RIN 5 , FOR Q 11 
DQDTN = 914.39 DMDT = 160.58 

RIN = 51.22 AIN = 113552.84 



VIN = 2576923.08 DVH20 = 35853817.74 

PDL = 31.04 PDLTOT= 1882.36 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

670.00 657.40 653.70 542.95 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.013799 0.004045 0.275825 0.003231 0.000000 

DQDTRO = 1.26 DQDTR3 = 512.88 DQDTN = 914.39 
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NUMBER OF ITERATIONS ARE-.RIN 5 ,FOR Q 10 
DQDTN = 267.75 DMDT = 47.02 

RIN = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 12.55 PDLTOT= 30.61 

THE TEMPERATURE PROFILE (T0.T2,T3 ,T4,T5) IS: 

590.00 584.37 583.72 540.87 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.019363 0.004312 0.374875 0.003231 0.000000 

DQDTRO = 2.92 DQDTR3 = 153.44 DQDTN = 267.75 

NUMBER OF ITERATIONS ARE:RIN 5 ',FOR Q 10 
DQDTN = 334,67 DMDT = 58.77 

RIN = 48.07 AIN = 107405.63 

VIN = 2307692.31 D’TH'^O = 453004.31 

PDL = 14.97 PEL ?T= . 39.38 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

600.00 593.95 592.52 541.08 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.018230 0.004276 0.356389 0.003231 0.000000 

DQDTRO = 2.74 DQDTR3 = 190.35 DQDTN = 334.67 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 10 
DQDTN = 405.52 DMDT = 71.21 

RIN = 48,53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 17.36 PDLTOT= 52.50 

THE TEMPERATURE PROFILE (TO ,T2,T3 ,T4 ,T5) IS: 

610.00 603.03 601.31 541.31 540.00 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.017301 0.004241 0,340714 0.003231 0.000000 

DQDTRO = 2.56 DQDTR3 = 229.42 DQDTN = 405.52 

NUMBER OF ITERATIONS ARE : RIN 5 , FOR Q 10 
DQDTN = 480.27 DMDT = 84.34 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.72 PDLTOT= ' , 76.82 

THE TEMPERATURE PROFILE (T0,T2,T3,T4.,T5) IS: 

620.00 612.11 610.09 541.55 540.00 

THE MATRIX VALUES A,B,C,D,E ARE : 

0.016516 0.004206 0.327048 0.003231 0.000000 

DQDTRO = 2.36 DQDTR3 = 270.72 DQDTN = 480.27 



NUMBER OF ITERATIONS ARE:RIN 5 .FOR Q 
DQDTN = 559.01 DMDT = 98.17 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.05 PDLTOT= 100.51 
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expansion of a fixed mass of Helium and the decrease in 

density of the air volume between the envelopes: 

mf(lbm) = q(BTU)/20500(BTU/lbm) 

change in lift (Ibf) = (To - Tinf) 130,000 

Tinf 

+ (W - Vjn) (T ..--. Tinf). 

Tinf 

'When steam injection is used to heat the inner volume, 
significant lift is produced by the volume the saturated 
steam displaces. Since the weight of the water is unchanged 
when it is vaporized, the change in lift produced by the 
steam is: 

change in lift (steam) (Ibf) = change in 
where is the increase in steam volume from initial 

(T^j^^) to heated (Tq) conditions. This volume is calculated 
for a specified temperature from: 

''steam * ''in^ R/ [ 18 . 02 (spv) p-RT) 
using the ideal gas equation for Helium and steam table 
values from Subroutine Volume for the steam. 

The loss of air spacing between the envelopes due to the 
additional inner volume from the steam is not accounted for 
by the program. 

Output quantities are filed for use by the IBM graphics 
program DISSPLA. 
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This concept has several significant disadvantages: 

- The double envelope will be heavier than conventional 
airship envelopes, although this factor is mitigated by 
not requiring additional ballonet envelopes. 



- A detailed design for the suspension and shape control 
of the envelopes has not been completed or tested. 

- Hot air ducting and/or a liquid fuel furnace will add 
■ weight to the basic airship. 

- A steam injection heating method requires a lot of water 
(14000 lb), energy to heat that water (7.32x10^ BTU for 
a 50 deg F temperature rise = 360 Ibm fuel) , pays a 
weight penalty for the liquid film on the inner envelope 
(about 7000 lb) , and requires some type of recovery 
system. 

- Availability, toughness, and cost of a low emissivity 
fabric coating are unknown factors. 

- Some time is required for heatup or cool down during 
takeoff and landing operations where load changes are 
significant, depending on heating method used and heat 
source . 

Overall, the promise of operational flexibility and 
capability should lead to an expanded study of this concept, 
with the effort and financing required to accomplish large 
scale testing to prove the estimated and unknown design 
factors, and solve those engineering and conceptual problems 
which will arise. 
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DQDT=DQDTN 

IF (UI.EQ.0.0) GO TO 500 
GO TO 600 
1000 CONTINUE 

C EQUIVALENT FUEL BURN RATE AND LIFT CHANGE VALUES ARE 
C DIRECTLY COMPUTED 

DMDT=DQDTN*3600 . 0/20500 . 0 

DLIFT=130000 .0*( (DELT/TI) )+( (TMEAN-TI)/TI)*(VTOT-VIN)*RHOMN 
+*32,174 

C VALUES BASED ON STEAM INJECTION HEATING ARE COMPUTED 
PDL=DLIFT*100 .0/130000 . 0 

CALL VOLUME (TI,SPVI) 

VH20I=VIN*TI*1545 . 0/ (18. 02*SPVI*14 .7*144 . 0-1545 . 0*TI ) 

CALL VOLUME (T0,SPV0) 

VH2OH=VIN*T0*1545 . 0/ (18 . 02*SPV0*14 .7*144 . 0-1 545 . 0*T0 ) 
DVH20=VH20H-VH20I 

DLH20=DVH20*RH0I*32 . 174- (VH20H*(TMEAN-TI ) *RHOMN*32 . 174) /TI 

DLT0T=DLIFT+DLH20 

PDLTOT=(DLTOT*100 .0)/130000 . 0 

C PROGRAM OUTPUTS ARE PRINTED; FILES FOR GRAPHICS USE ARE WRITTEN 

WRITE (*,1110) M,J 
WRITE (*,1120) DQDTN,DMDT 
WRITE (*,1130) RIN, AIN 
WRITE (*,1140) VIN, DVK20 
WRITE (*,1150) PDL, PDLTOT 

PRINT*, 'THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS:' 

WRITE (*,1100) T0,T2,T3,T4,T5 

PRINT*, 'THE MATRIX VALUES A,B,C,D,E ARE:' 

WRITE (*,1160) A,B,C,D,E 

WRITE (*,1180) DQDTRO, DQDTR3 , DQDTN 

WRITE (1,1100) UI,T0,DELT, DMDT , PDL , PDLTOT , E IN 
WRITE ( 2 , 1 1 00 ) UI , TO , DELT , DQDTN , PDL , PDLTOT 
WRITE (3, 1100) T0,T2,T3,T4,T5,TI,EIN 
WRITE (4,1100) TO , DELT , EIN , DQDTN , DMDT , PDL 
1100 FORI-IAT(F8.2,T10,F8.2,T20,F8.2,T30,F8.2,T40,F8.2,T50,F8.2,T60,F8.2) 

1110 FORMAT (IX, 'NUMBER OF ITERATIONS ARE;RIN' ,12, ' ,FOR Q ',12) 

1120 FORMAT (IX, 'DQDTN = ',F10.2,' DMDT = ',F10.2) 

1130 FORMAT (IX, 'RIN = ',F10.2,' AIN = ',F11.2) 

1140 F0P1IAT(1X, 'VIN = ',F10.2,' DVH20 = ',F11.2) 
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SECTION 1: VALUES OBTAINED FROM THE BASIC PROGRAM 



INPUT THE EMISSIVITIES OF THE INNER AND OUTER 
ENVELOPES, EIN AND EOUT 

. 1 , .9 

INPUT U INFINITY (UI,KTS), AND T INF (DEG R) 

50,540 

NUMBER OF ITERATIONS AREtRIN 5 ,FOR Q 8 
DQDTN = 41.16 DMDT = 7.23 

RIN = 45.75 AIN - 102801.46 

VIN = 2115384.62 DVH20 = 25472.64 

PDL = 2.57 PDLTOT= 3.98 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

550.00 548.80 548.61 540.29 540.16 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.031986 0.004468 0.549943 0.003231 0,003923 

DQDTRO = 3.54 DQDTR3 = 26.04 DQDTN = 41.16 



NUMBER OF ITERATIONS ARE: RIN 5 .FOR 0 9 



DQDTN = 89.71 DMDT = 

RIN = 46.22 AIN = 

VIN = 2153846.15 DVH20 = 

PDL = 5.14 PDLTOT= 

THE TEMPERATURE PROFILE (TO 
560'.00 557.76 557.36 

THE MATRIX VALUES A,B,C,D,E 
0.025960 0.004427 0.471630 

DQDTRO = 3.40 DQDTR3 = 



15.75 

103737.93 

79283.17 

9.50 

,T2,T3,T4,T5) IS: 

540.64 540.35 

ARE: 

0.003231 0.003923 

54.26 DQDTN = 89.71 



NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 9 
DQDTN = 143.07 DMDT = 25.12 

RIN = 46.68 AIN = 104666.37 

VIN = 2192307.69 DVH20 = 121494.44 

PDL = 7.67 PDLTOT= 14.32 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

570.00 566.80 566.17 ' 541.02 540.56 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.022918 0.004388 0.428758 0.003231 0.003923 

DQDTRO = 3.24 DQDTR3 = 84.42 DQDTN = 143,07 



NUMBER OF ITERATIONS ARE.-RIN 5 ,FOR Q 10 
DQDTN = 200.58 DMDT = 35.22 

RIN = 47.15 AIN = 105586.99 

VIN = 2230769.23 DVH20 = 192264.72 

PDL _= 10.18 PDLTOT= 20.63 

THE TEMPERATURE PROFILE (TO ,T2,T3,T4,T5) IS: 

580.00 575.87 575.00 541.44 540.79 
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APPENDIX B 



PROGRAM SAMPLE OUTPUTS 

The following program outputs are used to derive the 
values quoted in the results and conclusions chapters. 
These are divided into the following sections: 

- Section 1: Values obtained from the basic program. 



- Section 2: 


Values obtained for T 2 = Tq by setting "A" 
to zero in the program. 


- Section 3: 


Values obtained for T 2 = T^ by setting "C" 
to zero in the program. 


- Section 4 : 


Values obtained for T^ = by setting ”E" 

to zero in the program. 


- Section 5: 


Values obtained for increased envelope heat 
loss resistance by factors of 2 and 10. 



Units applicable to these data are: 

- Distances/lengths in ft. 

- Volumes in ft^ ; Areas in ft^ . 

- Heat flux in BTU/sec. 

- Fuel flow rate in Ibm/hr. 

- Temperature in degrees R. 
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1150 FORMATdX, 'PDL = ',F10.2,‘ PDLTOT= ',F11.2) 

1160 F0RMAT(1X,F3.6,T11,F8.6,T21,F8.6,T31,F3.6,T41,F8.6) 

1180 FORMAT(1X,'DQDTRO = ',F7.2,' DQDTR3 = ',F7.2,' DQDTN = ',F7.2/) 

C THE PROGRAM IS RUN FOR OTHER INPUT VALUES IF DESIRED 
1200 CONTINUE 
PRINT* 

PRINT* 

. PRINT*, 'DO YOU WANT TO ENTER OTHER VALUES OF UI,TI,EIN OR EOUT? ' 
PRINT*, '(1=YES, 2=N0)' 

READ*, Z 

IF (Z.EQ.1.0) GO TO 10 
STOP 
END 

SUBROUTINE VOLUIffi (T,SPV) 

THIS SUBROUTINE USES A SERIES OF LINEAR APPROXIMATIONS TO ' 
COMPUTE THE SPECIFIC VOLUME OF STEAM AT 1 ATMOSPHERE AT TEMP=T 

IF (T.LE. 530.0) THEN 

SPV=2445. 0+(T-500. 0)*(367. 7-2445. 0)/30 

ELSEIF (T.LE. 560.0) THEN 

SPV=867.7+(T-530.0)*(350.0-867.7)/30 

ELSEIF (T.LE. 590.0) THEN 

SPV=350. 0+(T-560. 0)*(157. 17-350. 0)/30 

ELSEIF (T.LE. 620.0) THEN 

SPV=157. 17+(T-590. 0)*(77. 23-157. 17)/30 

ELSEIF (T.LE. 650.0) THEN 

SPV=77. 23+(T-620. 0)*(40. 95-77. 23)/30 



ELSE 

SPV=40. 95+(T-650.0)*(23. 15-40. 95)/30 
ENDIF 
END 
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tabulated values at an intermediate temperature. A series 
of linear interpolations in Subroutine Volume is used to 
evaluate the specific volume of steam at a given temperature 
based on tabulated values from reference 12. Standard 
atmospheric pressure (14.7 psi) is used in all calculations. 

Since the linearized system of equations is of the 
simple form 
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a simple algebraic solution is used in the program for q, 
and then a substitution of new values is used to evaluate 
the temperature for the new iteration. 

Special computations are made if airship velocity or 
envelope emissivities are zero to avoid division by zero 
problems in the program. Results for airship velocities 
less than ten knots (but not zero) are less than accurate 
values due to the mixed convection conditions over the outer 
envelope . 

Equivalent fuel flow rate is based on 20500 BTU of 
useable heat energy per pound of fuel, which is obtainable 
from petroleum base fuels. Increase in lift due to heating 
without steam injection is based on the sum from the volume 
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THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004350 0.390713 0.003231 0.003923 

DQDTRO = 3.01 DQDTR3 = 131.32 DQDTN = 227.01 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 6 
DQDTN = 296.14 DMDT = 52.01 

RIN, = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 13.02 PDLTOT= 30.98 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

590.00 590.00 588.72 542.12 541.16 

THE MATRIX VALUES A,B,C,D.E ARE: 

0.000000 0.004312 0.368715 0.003231 0.003923 

DQDTRO = 2.82 DQDTR3 = 169.74 DQDTN = 296.14 



NUMBER OF ITERATIONS ARE:RIN 5 -FOR 0 6 



DQDTN = 369.62 DMDT = 

RIN = 48.07 AIN = 

VIN = 2307692.31 DVH20 = 

PDL .= 15.52 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

600.00 600.00 598.42 

theTmatrix values a,b,c,d,e 

0.000000 0.004276 0.350756 

DQDTRO = 2.62 DQDTR3 = 



. 64.91 
107405.63 
453004.31 
39.76 

,T2,T3,T4,T5) IS: 

542.64 541.45 

areT 

0.003231 0.003923 

210.61 DQDTN = 369.62 



NUMBER OF ITERATIONS ARE-.RIN 5 .FOR 0 6 



DQDTN = 447.46 DMDT = 

RIN = 48.53 AIN = 

VIN = 2346153.85 DVH20 = 

PDL = 17.98 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

610.00 610.00 608.10 
THE MATRIX VALUES A,B,C,D,E 
0.000000 0.004241 0.335507 

DQDTRO = 2.41 DQDTR3 = 



78.58 

108304.04 

655969.36 

52.85 

,T2,T3,T4,T5) IS: 

543.20 541.76 

ARE: 

0.003231 0.003923 

254.02 DQDTN = 447.46 



NUMBER OF ITERATIONS ARE:RIN'5 ,F0R Q 6 
DQDTN = 529.68 DMDT = 93.02 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 20.41 PDLT0T= 77.02 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

620.00 620.00 617.77 543.79 542.08 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.000000 0.004206 0.322198 0.003231 0.003923 

DQDTRO = 2.19 DQDTR3 = 300.06 DQDTN = 529.68 
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THE MATRIX VALUES A,B,C,D,E ARE: 

0.012216 0.004350 0.000000 0.003231 0.003923 

DQDTRO = 3.37 DQDTR3 = 0.00 DQDTN = 1688.06 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 8 
DQDTN = 2188.42 DMDT = 384.31 

RIN = 47.61 AIN = 106500.01 

VIN = 2269230.77 DVH20 = 333175.90 

PDL = 12.02 PDLTOT= 30.21 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS; 

590.00 565.09 555.66 555.66 548.58 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.011398 0.004312 0.000000 0.003231 0.003923 

DQDTRO = 3.27 DQDTR3 = 0.00 DQDTN = 2188.42 



NUMBER OF ITERATIONS ARE;RIN 5 .FOR 0 8 



DQDTN = 2705.01 DMDT = 

RIN = 48.07 AIN - 

VIN = 2307692.31 DVH20 = 

PDL = 14.42 PDLTOT= 

THE TEMPERATURE PROFILE (TO 

600.00 570.92 559.35 

THE MATRIX VALUES A,B,C,D,E 
0.010763 0.004276 0.000000 

DQDTRO = 3.17 DQDTR3 = 



475.03 

107405.63 

453004.31 

39.00 

,T2,T3,T4,T5) IS; 

559.35 550.61 

ARE: 

0.003231 0.003923 

0.00 DQDTN = 2705.01 



NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 8 
DQDTN = 3235.65 DMDT = 568.21 

RIN = 48.53 AIN = 108304.04 

VIN = 2346153.85 DVH20 = 655969.36 

PDL = 16.81 PDLTOT= 52.20 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

610.00 576.87 563.15 563.15 552.69 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.010249 0.004241 0.000000 0.003231 0.003923 

DQDTRO = 3.07 DQDTR3 = 0.00 DQDTN = 3235.65 

NUMBER OF ITERATIONS ARE:RIN'5 ,FOR Q 9 
DQDTN = 3778.79 DMDT = 663.59 

RIN = 48.98 AIN = 109195.42 

VIN = 2384615.38 DVH20 = 1071527.20 

PDL = 19.17 PDLTOT= 76.66 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

620.00 582.93 567.03 567.03 554.82 

THE MATRIX VALUES A,B,C,D,E ARE; 

0.009818 0.004206 0.000000 0.003231 0.003923 

DQDTRO = 2.95 DQDTR3 = 0.00 DQDTN = 3778.79 

NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 9 
DQDTN = 4333.20 DMDT = 760.95 

RIN = 49.44 AIN = 110079.94 
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The radiation equations are basically analytical, and 
should be accurate to within the uncertainties in 
emissivities and temperatures. 

■Overall unknown factors include additional conductive 
and convective losses in the lower portion of the airship 
due to the anchoring of the inner envelope, and losses in 
ducting and structural components. 

Although the combination of these uncertainty factors 

is : 

( 1 . 2 ) ( 1 . 4 ) ( 1 . 5 ) ( 1 . 25 ) = 3.15 

the major factor is the convection between the envelopes. 

It is believed that an airship constructed as conceived will 
have convective heat losses well within a factor of two 
times these results. Radiative heat losses, restricted to 
the assumptions stated, will be within a few percent of 
these results. 
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DELTA T (DEG F) 

Figure 5: Eq Fuel Flow vs Delta 



When this linear system is solved, these new solutions 
are resubstituted into the equations to find the heat 
transfer coefficients and radiative heat transfer values for 
the new iteration. This entire problem is iterated until 
convergence for q is obtained. 

Since several of the "constants" in the linearized 
system depend on the temperatures h^^, h^ , and k^) and 

particularly since the radiative heat transfer terms are 
fourth order in temperature, convergence of this iteration 
is expected to be a major problem. A close initial guess of 
the temperature values, however, produces convergence up to 
critical envelope emissivity values of more than 0.9. This 
is considered an adequate range for emissivities , and an 
adequate solution for this model. 

E. THE COMPUTER PROGRAM 

Program HOTSHIP utilizes these equations and completes 
the iterations with a convergence criterion of q within 0.1 
BTU/sec. From input values of airship velocity, Tq 

and envelope emissivities, the temperature profile and heat 
loss rate are computed. From these values, the increase in 
lift and equivalent fuel burn rate are directly computed. 

Properties of air and Helium with little temperature 
dependence are evaluated by a linear interpolation from 
tabulated values obtained from reference 8. The ideal gas 
equation is used in density change computations based on 
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NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 11 
DQDTN = 545.87 DMDT = 95.86 

RIN = 49.44 AIN = 110079.94 

VIN = 2423076.92 DVH20 = 1481890.78 

PDL = 22.24 PDLTOT= 100.51 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS ; 

630.00 621.33 619.06 543.91 542.14 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.015939 0.004172 0.316478 0.003231 0.003923 
DQDTRO = 2.16 DQDTR3 = 308.42 DQDTN = 545.87 

NUMBER OF ITERATIONS ARE-.RIN 4 ,FOR Q 11 
DQDTN = 626.30 DMDT = 109.98 

RIN = 49.89 AIN = 110957.85 

VIN = 2461538.46 DVH20 = 2272776.29 

PDL = 24.55 PDLTOT= 143.83 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4 ,T5) IS: 

640.00 630.42 627.83 544.48 542.46 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.015341 0.004139 0.305484 0.003231 0.003923 

DQDTRO = 1.95 DQDTR3 = 353.46 DQDTN = 626.30 

NUMBER OF ITERATIONS ARE:RIN 4 ,FOR Q 11 
DQDTN = 710.59 DMDT = 124.79 

RIN = 50.33 AIN = 111828.94 

VIN = 2500000.00 DVH20 = 4419910.85 

PDL = 26.83 PDLTOT= 257.35 

THE TEMPERATURE PROFILE (TO ,T2,T3 ,T4 ,T5) IS; 

650.00 639.50 636.58 545.08 542.79 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014808 0.004107 0.295_4_22 0.003231 0.003923 

DQDTRO = 1.73 DQDTR3 = 400.86 DQDTN = 710.59 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 11 
DQDTN = 798.84 DMDT = 140.28 

RIN = 50.78 AIN = 112694.08 

VIN = 2538461.54 DVH20 = 8108853.38 

PDL = 29.08 PDLTOT= . 449.13 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS: 

660.00 648.58 645.32 545.72 543.13 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.014327 0.004075 0.286114 0.003231 0.003923 

DQDTRO = 1.50 DQDTR3 = 450.71 DQDTN = 798.84 

NUMBER OF ITERATIONS ARE; RIN 5 ,FOR Q 11 
DQDTN = 891.11 DMDT = 156.49 

RIN = 51.22 AIN = 113552.84 

VIN = 2576923.08 DVH20 = 35853817.74 

PDL = 31.30 PDLTOT= 1876.12 

THE TEMPERATURE PROFILE (T0,T2,T3,T4,T5) IS; 

670.00 657.64 654.04 546.38 543.50 
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NUMBER OF ITERATIONS ARE;RIN 5 ,FOR Q 12 
DQDTN = 598.22 DMDT = 105.05 

RIN = 51.22 AIN = 113552.84 

VIN '= 2576923.08 DVH20 = 35853817.74 

PDL = 31.43 PDLTOT= 1872.96 

THE TEMPERATURE PROFILE (T0,T2,T3 ,T4,T5) IS: 

670.00 660.84 636.64 561.68 542.35 

THE MATRIX VALUES A,B,C,D,E ARE: 

0.015348 0.040446 0.298304 0.032314 0.003923 

DODTRO = 1.17 DQDTR3 = 346.92 DQDTN = 598.22 

NUMBER OF ITERATIONS ARE:RIN 5 ,FOR Q 13 
DQDTN = 657.62 DMDT = 115.48 

RIN = 51.66 AIN = 114405.53 

VIN = 2615384.62 DVH20 = *********** 

PDL = 33.62 PDLTOT= • -807.37 

THE TEMPERATURE PROFILE (TO ,T2 ,T3 ,T4,T5) IS: 

680.00 670.19 643.79 563.83 542.58 

THE MATRIX VALUES A,B,C,D,E ARE : 

0.014937 0.040145 0.289940 0.032314 0.003923 

DQDTRO = 0.91 DQDTR3 = 381.83 DQDTN = 657.62 
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APPENDIX C 



COMPUTATION OF AVAILABLE ENGINE HEAT 

From publications on powerplant characteristics, such as 
the engine manufacturer's specifications or Aviation Week , 
some basic trends for turboshaft and turboprop engines such 
as would be used for propulsion on this airship concept are: 
exhaust gas temperature = 600 deg C = 1100 deg F 
mass flow rate (Ibm/sec) = shp/110 (Ibm/sec) 

It is estimated that two 1000 shp engines (at least) 
would be used on the concept design. If the exhaust 
temperature could be reduced by 500 deg F in a heat exchange 
system without significantly effecting engine performance, 
then the exhaust heat available from each engine is: 
q = (mass flow rate) Cp^j^j. (500) 

= (1000/110) (7.7/32.174) (500) BTU/sec 
= 1088 BTU/sec 

Certainly almost twice as much heat (a 1000 deg F drop 
in temperature) is available if the exhaust is ducted 
directly into an air volume. 

Some additional heat is available from cooling the 
engine casing, but an estimate of quantity is not completed. 

Approximately one third of these heat values should be 
available at low or idle power, and about one half is 
expected to be available at moderate power. 



81 



APPENDIX D 



CALCULATION OF VALUES FOR STEAM INJECTION 

1. Estimation of Water Film Weight 

From Gebhart (reference 10, pp. 435-436) the average 
thickness of a liquid water film on a large vertical plate 
will be about 0.005 inches. Based on this thickness, 
Weightj^20 film “ Area (thickness) ^ g 

= 1.5xl0^ft^ (.005 ft/12) (60.81b/ft^) 

= 3800 Ibf 

From the additional film thickness at the top and bottom 
of the envelope, the total water film weight is expected to 
double this value, or as much as 7000 Ibf. 

2 . Mass of Steam Evaporated 

At 110 deg F, the mass of water vapor in 2xl0®ft^ of 
saturated Helium is: 

™H20 ^ ^ 2x10® /265.1 = 7544 Ibm 

3 . Energy to Evaporate Steam 

For the mass of water evaporated at 110 deg F, 
heat required = (970BTU/lbm) (75441bm) 

= 7.32 X 10® BTU 

or in terms of equivalent fuel (kerosene) , 

^kerosene "" 7.3 2x10® BTU/ (20500 BTU/lbm) 

= 360 Ibm 

If engine heat of 2000 BTU/sec is used, then the time 
required to evaporate this much steam is: 

time = [7.32xl0®/2000] (1 hr/3600 sec) = 1 hr 
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